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ABSTRACT 

We analyze the effect of internal degrees of freedom on the movability prop- 
erties of localized excitations on nonlinear Hamiltonian lattices by means 
of properties of a local phase space which is at least of dimension six. We 
formulate generic properties of a movability separatrix in this local phase 
space. We prove that due to the presence of internal degrees of freedom of 
the localized excitation it is generically impossible to define a Peierls-Nabarro 
potential in order to describe the motion of the excitation through the lattice. 
The results are verified analytically and numerically for Fermi-Pasta-Ulam 
chains. 
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Recently localized breatherlike excitations were discovered to exist in sev- 
eral different Hamiltonian lattices in one and two dimensions [0,0,0,0], 
They are self- localized (no disorder) and appear in nonlinear lattices 
- thus we name them nonlinear localized excitations (NLEs). For certain 
systems it was possible to create moving NLEs Consequently the idea 

arose to describe their motion in a Peierls-Nabarro- Potential (PNP) pfl , [|1C[| . 



1J1,[12[ related to the PNP of kinks |[L3|,|14||. Numerical simulations strongly 



support the existence of a PNP-related phenomenon in Fermi-Pasta-Ulam 
systems and Klein-Gordon systems fl6| . However as we show below it is 
generically impossible to define a PNP for NLEs. 

The NLE solutions are nontopological, i.e. no special structure of the 
underlying many-particle potential is required. The only condition is to have 
nonlinear terms in the potential. One can perform stability analysis and show 
that if the NLE is localized enough (in practice it will contain only a few 
particles which are involved in the motion) then generically all Hamiltonian 
lattices will exhibit families of stable time-periodic NLE solutions Hi , fl7l , [[IS 



Hereafter we will call these stable periodic NLEs elliptic NLEs to emphasize 
their stability property (in a Poincare mapping they would appear as stable 
elliptic fixed points |J). One can view the NLE as a solution of a reduced 
problem where only M particles are involved in the motion, the rest of the 
lattice members are held at their groundstate positions. We showed that 
many frequency NLEs can be excited by perturbing the elliptic NLEs and 
that thus NLE solutions are motions on M-dimensional tori in the phase space 
of the reduced problem and in the corresponding local subspace of the phase 



space of the full system p|, p7[| . Besides these stable NLE solutions unstable 
periodic NLEs exist. Their feature is that certain local perturbations destroy 
the unstable NLE or cause it to move [p|| , P~5[| , [1T9B . Hereafter we will call 
them hyperbolic NLEs. If one calculates the energy density distribution e/ 
for the NLE solutions, one can define the position of the energy center of 
the distribution by X E — • e i/(J2i e i)- F° r a given system the elliptic 
NLE solution yields X E = l (i.e. centered on a lattice site Z ) and the 
hyperbolic NLE solution X E = h + 0.5 (i.e. centered between lattice sites 
h and l\ + 1) or vice versa || , , |19|| . Here l,lo,h denotes lattice sites 



and the lattice spacing is 1. Both elliptic and hyperbolic NLEs as well as 
certain stable subclasses of their perturbations obey a symmetry during their 
whole evolution namely that the evolution of the NLE part for x < X E is 
symmetric (or antisymmetric) to the evolution of the NLE part for x > X E . 
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This symmetry is just the manifestation of Hamiltonian character of motion 
combined with the discrete translational symmetry of the lattice. 

The writing down of a certain PNP for the collective coordinate which 
describes the motion of the NLE is conceptually equivalent to the problem 
of a pendulum. The PNP-barrier Ap^r is then intimately connected with the 
energy that is required to overcome the separatrix of the pendulum. This 
separatrix separates oscillating pendulum solutions from rotating ones [pCfl . 
The PNP frequency ujpn is essentially the pendulum frequency for infinitely 
small amplitudes. 

To describe a periodic elliptic NLE we need to introduce one degree of 
freedom. We will work in the action-angle phase (J, 9) space and name 
this degree of freedom J\. Its corresponding frequency will be uj\ = Q\ = 
dH/dJ\. Here H denotes the full Hamiltonian of the lattice. We assume that 
there exists a certain transformation between the original variables (positions, 
momenta) and the actions and angles. This does not imply integrability of 
the system as well as it does not imply the inverse. Since our NLE solutions 
are regular solutions (at least on moderate time scales) there is no need in 



introducing stochasticity (cf. |T7j for details). Because of the symmetry of 
the elliptic NLE the NLE will be stationary (nonmoving) for any value of J\ in 
the whole range of its existence. To excite a moving NLE we have to excite an 
additional degree of freedom J3. Exciting J3 we destroy the symmetry of the 
elliptic NLE. But since it is always possible to perturb the NLE conserving 
the symmetry, we have to include an additional symmetryconserving degree 
of freedom J 2 into the consideration. Thus we end up with the simplest 
generic case of a Hamiltonian problem with three degrees of freedom: 

OH 

H = H(J 1 ;J 2 ;J 3 ) , = 0, = — , i = 1, 2, 3 . (1) 

According to our notation i = 3 labels the symmetrybreaking degree of 
freedom. If it is excited strongly enough we expect to hit a separatrix which 
separates stationary NLEs from moving ones. We will name this separatrix 
movability separatrix. All three degrees of freedom are of local character, 
especially they can be well defined in the reduced problem for the NLE. 
Since we can consider the NLE excitation at (or between) any lattice site(s) 
we thus study the local character of a movability separatrix which is also 
defined for the infinite system. The movability separatrix for the full system 
is just a periodic continuation of the local movability separatrix. 
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Let us state the general condition for the movability separatrix we are 
looking for. Since on the movability separatrix a trajectory will for inifmite 
times asymptotically reach a hyperbolic state (which is nothing else then 
the hyperbolic NLE and its symmetric perturbations) the corresponding fre- 
quency of the 3d degree of freedom 

U3 = W = /(Jli Ja; Js) (2) 

has to vanish on the movability separatrix i.e., 

/(Ji;J 2 ;J 3 ) = , (3) 

which implies an equation for a surface in the three-dimensional subspace of 
the actions (Ji; J 2 ; J3). We can always eliminate J 2 using the expression for 
the energy E = H(J\\ J 2 ; J3), so that |] yields: 

f(Ji,J 2 ;J 3 )=f(E;J 1 ;J 3 ) = . (4) 

From [| it follows that there exists a critical value for J 3 on the movability 
separatrix: 

J s 3 =g(J 1 ;J 2 ) = ~g(E-J l ) • (5) 

The only possibility of introducing the PNP would be to use the relation 
between the potential of a pendulum and its critical value for the action as 
well as the frequency of small amplitude oscillations: 

to PN = f(J i; J 2 ;J 3 = 0) = f(E;J i; J 3 = 0) . (6) 

It is very important to note that if / from ^ or / from ^ depend on ( Ji, J 2 ) 
or (E; J\j respectively then the same fact holds for ujpn in |^ as well as for 
J| in ^. As we immediately recognize a PNP would be different for different 
(E; Ji) because of the generic dependence of the PNP parameters on the 
values of E and J\ in It is not only that we would obtain different PNPs 
by varying the energy. Even for a fixed energy different PNPs occur because 
of the dependence of the right-hand sights in on J x . 

Let us discuss special nongeneric cases: i) The Hamiltonian separates in 
the actions in the following way: 

H{Ju J 2 ; Js) = H 12 {J X - J 2 ) + ^3(^3) (7) 
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Then according to its definition u)% depends only on J3. Thus the value 
of Jf becomes independent on (Ji; J 2 ) or (E; J\) and a unique PNP can 
be immediately associated with the term H 3 (J 3 ) in [7|. ii) A more subtle 
nongeneric case appears if no separation holds but the frequency u 3 is only 
a function of energy E. In this PNP can be introduced which would 

depend on the energy of the NLE. 

Let us apply the results from above to a class of systems where moving 
NLEs were detected J7|,| 
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H 



£ ( -Pf + V{X l - x w 



V(x) 



l -Cx* + \x* 



(8) 
(9) 



These systems belong to the class of Fermi- Past a- Ulam systems pl |. Pi and 
Xi are momentum and position of the Z-th particle respectively. The param- 
eter C regulates the strength of the quadratic terms. For C — > 00, E = const 
Eq.^| becomes the well-known linear atomic chain, which is integrable and has 
no NLE solutions. For C — > 0, E = const Eq|| becomes a highly nonlinear 
nonintegrable atomic chain. All properties of |8] can be obtained by fixing the 
energy e.g. at E — 1 and varying C. All solutions for other energies can be 
obtained by proper scaling of the times, displacements and the parameter C: 
if {Xi(t; E = 1; C)} is a solution of |], then {Xi(t; E; C} is a solution for the 
energy E = A -4 , parameter C = C/\ 2 and Xi(t) = A _1 X(A _1 /f:). Let us first 
discuss the case C = 0. Then an even simpler scaling holds - it is enough to 
study the system at one given energy e.g. E = 1 and through the above de- 
scribed scaling all solutions for other energies are obtained. The elliptic NLE 
solution is the well known even parity mode P|,fll5]l. It is centered between 
two particle sites (X E = I + 0.5) and four particles are essentially involved 
in the motion. Its amplitude distribution can be qualitatively indicated by 



(■ " T 



i • •)■ More precisely the scaled absolute values of the amplitudes 
in decreasing order read: 1, 0.16579, 0.00048, ... . No exact compacton 
structure is observed as it was wrongly claimed in ]TIJ because of a calcula- 

The frequency of the elliptic NLE for E 



tion error in eq.13 of JTT 
ujAE 



1 is 



1) = 1.760 ± 0.0018. The hyperbolic NLE solution is known as the 
odd parity mode |I5||,|119|. It is centered on a particle (X E = I) and essen- 
tially three particles are involved in the motion. Its amplitude distribution is 



5 



roughly (• • j J, • •). More precisely the scaled absolute values in decreasing 
order read: 1, 0.52304, 0.02305, ... . The frequency of the hyperbolic NLE 
is found to be u h (E = 1) = 1.751 ± 0.0018. 

Let us mention an important property of || Besides the energy conser- 
vation law this system conserves the total mechanical momentum: J2i Pi — 
const. It is sufficient to study the system in the center of mass frame, so 
that the total momentum vanishes and the center of mass doesn't move. All 
other cases can be obtained by a Galilean boost in §. Since the NLE solution 
is localized, the total mechanical momentum outside the NLE is zero. Thus 
it has to be zero inside too, i.e. our NLE solutions have to obey mechanical 
momentum conservation, at least approximately. The consequence is that 
the elliptic NLE (four particles) is described by 4 — 1 =3 degrees of freedom. 
That is exactly our simplest generic problem as described above. 

The properties of the perturbed elliptic NLE can be studied with Poincare 
mappings for symmetry-preserving perturbations, i.e. for J 3 = 0. Then we 
can consider a reduced problem where the particles outside the NLE are 
fixed at position zero. This fixed boundary doesn't break the momentum 
conservation because of the antisymmetry of the perturbed NLE. The result 
for the Poincare map is shown in Fig.l. The point in the middle of the map 
corresponds to the elliptic fixed point solution. All torus intersections inside 
the diamond-like structured torus correspond to stable two-frequency NLEs 
in the full system (1000 particles). Every torus in Fig.l corresponds to a 
certain triple of (J±] J<i\ Jj, = 0). The fixed point (elliptic NLE) is defined 
by (Ji; J2 = 0; J3 = 0). Thus we first arrive at the unambiguous result that 
a perturbation of the elliptic NLE preserving the symmetry leads to two- 
frequency NLE solutions (Ji; J 2 ; J3 = 0). This is similar to NLE properties 
in Klein-Gordon lattices P|,[^7|. 

Now we excite the third degree of freedom J3 7^ which destroys the 
symmetry of the elliptic NLE. We choose a path in phase space where P\(t = 
0) = -P Q (t = 0), P 2 (t = 0) = -P^{t = 0) = s, X x {t = 0) = -X_x(t = 0) = 
a and all other displacements/momenta are equal to zero at t — 0. The total 
energy is still E — 1. We work with 1000 particles. Here the elliptic NLE is 
chosen to be centered between the lattice sites / = and I = 1 respectively. 
The actions are some functions of the choosen path: J x = Ji(E;s;a) and 
^3 = J?,{.E] s; a). Especially we know that Js(E, s, a = 0) = 0. By increasing 
a we measure the time dependence of the energy center X^(t). The energy 
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density is defined by 

ei = \pf+ l -{v{x l -x l ^ 1 )+v{x l+l -x l )) . (10) 

Since X E {t) is independent of time for a = 0, we can hope that the en- 
ergy center will essentially couple only to ( J 3 ; 6 3 ) so that we can measure 
the frequency uj 3 . Indeed for a ^ X E {t) oscillates around its mean value 
of 0.5. There are modulations of this oscillation with the frequency Ui, 
but their amplitude is small and we clearly observe the frequency uj 3 = 
f(E; J\(E; s; a); J${E\ s; a)). For small values of a (< 10~ 4 ) the value of uj 3 
becomes independent of a, thus we can measure f(E; Ji(E; s; a = 0); J 3 = 0) 
which is nothing else than ujpn (cf. []). Especially for the elliptic NLE 
solution we find uo 3 = 0.343 ± 0.006. For two other tori within the dia- 
mond like torus in Fig.l (s = 0.1, s = 0.2) we find uo 3 = 0.391 ± 0.005 and 
uj 3 = 0.322 ± 0.005. Thus we find variations of uj 3 = f(E; J\\ J 3 = 0) for a 
fixed value of the energy (by varying J\) of at least 21%. Now let us increase 
a for a given value of s and monitor the time evolution of X^(t). It is shown 
in Fig. 2. In agreement with our expectations we find that with increasing 
a the frequency uj 3 decreases and the amplitude of the oscillations of X E (t) 
increases. At a threshold value of a = a s we clearly observe the crossing of 
the movability separatrix - the NLE escapes from its mean position. 

The properties of the movability separatrix are easy constructed. Because 
of the scaling property of the Hamiltonian for C = we find all solutions at 
other energies by proper scaling. Since the frequencies scale too, we imme- 
diately find that f(E; J\\ J 3 = 0) depends on the energy. Because we found 
strong (20%) variation of this frequency on the energy hypersurface (for con- 
stant energy) the Ji-dependence is also significant. Having f(E; J\] J 3 = 0) 
to be strongly dependent on E and J\ we find using that the same holds 
for the critical value of J 3 on the movability separatrix. Thus we see that 
our example is a generic case, and a PNP can't be constructed. 

If one considers C / (here C = 0.3) one rediscovers all the above 
statements. There are only quantitative changes - the dependence of u 3 on 
E and J\ becomes weaker. For large enough values of C (fixing the total 
energy) the frequencies UJ2 and u 3 can become resonant with the phonon 
band (which still doesn't prevent us from studying the movability separatrix 
on short time scales). For too large values of C the frequency w\ becomes 
resonant with the phonon band and the whole NLE solution then quickly 
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dissappears p2| , f23| , P , fT7 

Let us make some final comments. First our results demonstrate a clear 
way of studying and characterizing the movability properties of NLEs in 
terms of a movability separatrix in phase space. Secondly we find that 
generically no simple PNP can be introduced. The reason is the intimate 
connection between the 'translational' degree of freedom ( J 3 ) and the 'inter- 
nal' degrees of freedom (J 1; J 2 ) through the Hamilton function. Consequently 
the necessary energy supply to an elliptic NLE in order to cross the movabil- 
ity separatrix at a certain orbit can be positive, zero or negative depending 
on the choosen orbit on the movability separatrix. That is the reason why 
intuitive approaches to derive PNPs are sometimes even selfcontradictory: in 
r^] under assumption of separability property (our equation |?p of a discrete 



nonlinear Schrodinger equation a PNP is derived which is energy dependent, 
but that implies the nonseparability of the Hamiltonian. The results in the 



present paper disprove the conjecture in jnj, where it is predicted that for 
our equations |||] and C = freely moving NLEs exist , i.e. no PNP (no 
separatrix) should exist. Our Fig. 2 shows that the separatrix exists. A very 



interesting perturbation analysis was carried out in [10|],[|Ij] for a weakly 
perturbed integrable Ablowitz-Ladik lattice. The authors were able to show 
analytically that the NLE solution is described by the evolution of three in- 
ternal degrees of freedom, so that the movability separatrix can be analyzed 
analytically in their case. Finally we mention the treatment of a discrete 



sine-Gordon breather with a collective coordinate method p4| . There it was 
shown how to treat unambiguously a NLE with two degrees of freedom. Al- 
ready there it is clear that no unique PNP exists, i.e. the amplitude of the 
PNP is a function of E and the NLE amplitude. 

This work was supported in part (S.F.) by the Deutsche Forschungsge- 
meinschaft (F1200/1-1). 
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FIGURE CAPTIONS 



FIG.l 



Poincare intersection between the trajectory 



of system ||| for a reduced problem with fixed boundaries: 
C = 0, E = 1, P 1 = -P , X x = -X , P 2 = X 2 = -X. 

all other lattice members are fixed at position zero. 



FIG. 2 Time dependence of the center of energy of an NLE 
for different assymetric perturbations a (see text): 
solid line - a = 0.02; dotted line - a = 0.06; dashed line - a = 
long dashed line - a = 0.112; dashed-dotted line - a = 0.113. 



